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We present an update on our on-going project to compute hadronic observables for Nf = 2 flavours 
of O(a) improved Wilson fermions at small lattice spacings. The procedure to determine the lat- 
tice scale via the mass of the Omega baryon is described. Furthermore we present preliminary 
results for the pion form factor computed using twisted boundary conditions, and report on the 
implementation of a novel approach to determine the contribution of the hadronic vacuum polar- 
isation to the anomalous magnetic moment of the muon. 
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1. Introduction 

Studies of structural properties of mesons and baryons, encoded in their electromagnetic form 
factors, are the subject of many recent lattice calculations. However, systematic uncertainties aris- 
ing from a variety of sources appear to be much larger than for simpler quantities such as hadron 
masses and decay constants. As part of the CLS project 1 a set of ensembles is being generated, 
which will allow for a thorough investigation of systematic effects, such as lattice artefacts, finite- 
volume effects, and extrapolations in the light quark masses. The CLS project is based on proven 
and conceptually simple technology: simulations are performed using non-perturbatively 0(a) im- 
proved Wilson quarks and the Wilson plaquette action. In order to preserve the local structure of the 
lattice action, the link variables remain unsmeared. We employ the deflation-accelerated DD-HMC 
algorithm [1] for Nf = 2 flavours of light quarks. 

A potentially severe problem for any lattice simulation near the continuum limit is the ob- 
served sharp rise in the autocorrelation time of the topological charge [2]. This important issue 
has been addressed in several papers [3,4] and contributions to this conference [5]. In spite of the 
progress made in identifying the causes of this phenomenon, a rigorous and effective solution is still 
elusive. Therefore we have restricted the calculation of observables to a range of bare couplings 
for which the problem of very long autocorrelation times in the topological charge is not observed. 
A compilation of these ensembles is shown in Table 1. At our smallest value of the lattice spacing, 
a Ri 0.05 fm (i.e. at /3 = 5.5), where the problem is expected to be most severe, we have checked 
explicitly that the topological charge fluctuates around zero at a sizeable rate and produces a rea- 
sonably symmetric distribution [6] . We can therefore take confidence that the composition of our 
ensembles is not strongly biased and that our statistical errors are reliable. 

In this note we report on preliminary results for the pion electromagnetic form factor. 2 Another 
issue which has received a lot of attention recently, is the hadronic vacuum polarisation contribution 
to the muon's anomalous magnetic moment, a^ d . Lattice calculations of this quantity involve the 
determination of the vacuum polarisation amplitude, which, like hadronic form factors, depends on 
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Table 1: Simulation parameters and approximate values for the lattice scale and pion masses for those 
ensembles which show an acceptable tunnelling rate of the topological charge. The preliminary results 
presented in theses proceedings are based on the ensembles labelled "N", "E" and "F6". 
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2 Form factor calculations for baryons are presented in another contribution to this conference [7]. 
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Figure 1: Left: the effective mass in the £2-channel for the ensemble N5 (m K w 410MeV); Right: Sea 
quark mass dependence of mq for datasets N3 — N5. Red and black symbols denote the results from corre- 
lated and uncorrected fits, respectively. 

a momentum variable. Therefore, we discuss both the calculation of mesonic form factors and the 
determination of aj^ ad in these proceedings. 

2. Setting the scale 

Studying the scaling properties of observables and expressing dimensionful quantities in phys- 
ical units requires the computation of a reference quantity which sets the scale. The mass of the 
Q. baryon has emerged as a good candidate, since the Q. is stable in QCD, and its valence sector 
consists entirely of strange quarks. Moreover, ChPT studies [8] suggest a simple functional form 
for chiral extrapolations in the sea quark mass. 

In order to determine the masses of pseudoscalar and vector mesons, as well as octet and 
decuplet baryons, we used Jacobi-smeared sources [9], supplemented by HYP-smeared [10] link 
variables. In the pseudoscalar channel long and stable plateaus were observed, which allowed for 
a reliable determination of the mass via single-cosh fits to the correlation functions. In the vector 
channel and for baryons, the unambiguous identification of the ground state turned out to be more 
difficult. We therefore employed the procedure of [11], which is based on an ansatz that includes 
contributions from the first excited state. In Fig. 1 we show a typical effective mass plot in the 
£2-channel together with the fit result. 

In order to determine an interpolation in the valence quark mass must be performed. This 
requires knowledge of the hopping parameter, K s , which corresponds to the bare strange quark 
mass. We have determined K s using the procedure described in [6, 12]: Denoting the generic non- 
degenerate pseudoscalar and vector mesons by tkk and m* K , respectively, we interpolate the mass 
ratio (mK/m* K ) 2 as a function of (amjf ) 2 to the experimental value mK/m* K = 0.554. This yields the 
value of K s at each value of the sea quark mass. In a second step one interpolates aniQ in the valence 
quark mass to the value of K s . The resulting estimates of aniQ as a function of 1/ fc sea are shown for 
j3 = 5.5 in the right panel of Fig. 1. A mild dependence on the sea quark mass is observed. At this 
stage we have not yet attempted a chiral extrapolation in the sea quark mass to the physical value 
of m % . In order to set the scale at jS = 5.5 we therefore take the result for aniQ at the smallest sea 
quark, which yields = 0.053(1 )fm. 
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Our production runs for baryonic two-point functions on the other ensembles are not yet com- 
pleted. In order to compare form factors computed at different lattice spacings consistently, we 
employ a simple scaling relation, i.e. 

, / «ref 1/3=5.3 \ , 

a n|j3=5.3 = — -j a ^l/3=5.5 • C 2 - 1 ) 

V a ref|/3=5.5/ 

Here the scale a re f is defined via the kaon mass amg, evaluated at the (unphysical) reference point 
where m K jmx = 0.85 [12], which can be reached without performing a chiral extrapolation. At 
J3 = 5.3 this yields a ref = 0.0784(10) fm [12], while at j3 = 5.5 one obtains a ref = 0.0603( 15) fm [6]. 
Our result for at j3 = 5.5 and the scaling relation eq. (2.1) imply 

^1/3=5.5 = 0.053(l)fm, a a \ p=5 3 = 0.069 (2) fm. (2.2) 

At P = 5.5 the value of a& is smaller by more than 10% compared to a re f, which we used previously 
to convert our pion masses on the "N"-lattices into physical units (see Table 2 of ref. [6]). Using 
the value in eq. (2.2) the pion masses change to m % = 600, 510 and 410MeV on N3, N4 and N5, 
respectively. 



3. The pion electromagnetic form factor 



The electromagnetic form factor, defined by 

{^ + {Pf)\luy^u-\dy^d\% + {pi)) = (pf + p^Mq 2 ) 



(3.1) 



where q = pf — pt is the momentum transfer, encodes the distribution of electric charge inside the 
pion. Of particular interest is the charge radius, which is derived from the pion form factor at 
vanishing momentum transfer, i.e. 



(3.2) 



Lattice calculations of mesonic matrix elements are technically simpler than the corresponding 
quantities for the nucleon. Furthermore, the pion electromagnetic form factor receives no contri- 
butions from quark-disconnected diagrams. This opens the possibility to perform a precision test 
of lattice QCD, by comparing lattice estimates for (r 2 ) to the experimentally determined value. 
However, owing to the finite spatial volume the accessible range of momentum transfers Q 2 = —q 2 
is severely constrained, which presents a major obstacle for precise lattice determinations of (r 2 ). 
We have therefore employed flavour-twisted boundary conditions [14] in the valence sector, such 
that the momentum transfer satisfies [15] 



-Q z = q l = (pf-piY = [E K (pf)-E x (pi) 



p f + 



e, 



(3.3) 



We tuned the twist angles , 0/ so as to achieve a particularly fine momentum resolution near Q 2 = 
0. In order to improve the statistical signal we used stochastic Z2 x Z2 sources in the calculation 
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Figure 2: Left: pion form factor computed for a range of pion masses compared to the results of [16, 17]; 
Right: the charge radius as a function of the pion mass, compared with the results of the 2 + 1 -flavour 
simulation from [16]. 



of quark propagators. As in ref. [16] the pion form factor was extracted from a suitable ratio of 
correlators in which the renormalisation factor of the local vector current, Zy, drops out. 

In Fig. 2 we show our results for the pion electromagnetic form factor as a function of the 
squared Euclidean momentum transfer for the data sets N3, N4, N5 and F6. The corresponding 
pion masses in MeV are 600, 510, 410 and 290, respectively. Twisted boundary conditions were 
also employed in refs. [16, 17], which are shown together with our data. The much larger number 
of data points near Q 2 = allows us to perform a very accurate and largely model-independent 
determination of the pion charge radius, by computing the linear slope of f K (Q 2 ) in the vicinity 
of Q 2 = 0. In this way our estimate of the pion charge radius does not rely on a particular ansatz, 
such as vector dominance. As can be seen from the right panel in Fig. 2 we obtain statistically 
very precise results for (r~), in particular on the ensembles at j8 = 5.5. Statistics at our lightest 
quark mass (/3 = 5.3) will be increased. Our preliminary results for (r 2 ) are consistent with the 
determination of ref. [16], which uses Nf = 2+ 1 flavours of domain wall quarks, and the trend in 
our lattice data compares favourably with the experimentally determined value. 



4. Hadronic vacuum polarisation contribution to (g — 2)^ 

The anomalous magnetic moment of the muon, = ^ (g — 2)^, is among the most precisely 
measured quantities. Assuming the validity of the Standard Model, the experimental value for 
differs from the theoretical prediction at the level of 3.2 standard deviations [18]. An important 
ingredient is the leading hadronic contribution to vacuum polarisation, . This quantity is nor- 
mally determined via a semi-phenomenological approach based on the evaluation of a dispersion 
integral containing experimentally measured hadronic cross sections. Given the importance of 
for new physics searches, a first-principles determination of a|f d is highly desirable. Following 
refs. [19,20], a^ d can be computed on the lattice via the convolution integral 

< d = 4tt 2 ) 2 [ d<2 2 f(Q 2 ) (n(G 2 ) - 11(0) } , (4. 1) 
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Figure 3: Left: the vacuum polarisation amplitude computed at j3 = 5.3, L ~ 2.2fm, m n m 550MeV, using 
twisted boundary conditions. The solid curve is a Pade fit to the data; Right: the integrand of the convolution 
integral on the "F6" lattice. The area underneath the solid curve yields the value of a^ ad . 

where H(Q 2 ) is related to the vacuum polarisation tensor ILjy, i.e. 

n MV (2) = / d 4 W^') </ M (x)7 v (j)> = (q^qy-8^q 2 )U(q 2 ). (4.2) 

Here denotes the flavour-singlet vector current. A well-known problem is that the convolution 
function f{Q 2 ) is strongly peaked for momenta near m», which is an order of magnitude smaller 
than what can conventionally be realised on current lattice sizes. Twisted boundary conditions 
cannot be used in a straightforward manner, since the current-current correlator contains quark- 
disconnected diagrams, for which the effect of the twist cancels. 

To tackle this problem we have implemented the strategy outlined in [21]: The key observation 
is that the continuum limit exists for the individual quark-disconnected and connected contributions 
to IT. Moreover, the connected part can be re-interpreted in terms of flavour-non-singlet correla- 
tors, which are easily computed using twisted boundary conditions in the standard fashion. ChPT 
at NLO also predicts that the contribution from disconnected diagrams is suppressed by a factor 
— 1/10 relative to the connected ones [21]. 

Preliminary results for the vacuum polarisation amplitude extracted from the connected cor- 
relator in the two-flavour case are shown in Fig. 3. Twisted boundary conditions produce a much 
larger density of Q 2 values, which greatly stabilises the extrapolation to Q 2 = and the extraction 
of n(0), which enters eq. (4.1). Furthermore, owing to the larger reach in Q 2 the integrand of the 
convolution integral is much more tightly constrained (c.f. Fig. 3, right panel). Our lattice estimates 
for a]jf computed for Nf = 2 and also in the theory containing a quenched strange quark will appear 
elsewhere [22]. Eventually we will also include quark-disconnected diagrams, computed for the 
usual Fourier-momenta only, in order to verify that their contribution is indeed suppressed. 

5. Summary and outlook 

The ensembles generated as part of the CLS project allow for a comprehensive investigation 
of systematic effects for a variety of hadronic observables. Our fine lattice resolution will enable 
reliable continuum extrapolations and will prove useful to explore hadronic form factors at large 
momentum transfers. Twisted boundary conditions are indispensable for an accurate determination 
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of the pion form factor. Their use can also be extended to lattice determinations of ah , for which 
a straightforward application of this technique is not obvious. 

Acknowledgments: We thank our colleagues within the CLS project for sharing gauge ensembles. 
Calculations of correlation functions were performed on the dedicated QCD platform "Wilson" at 
the Institute for Nuclear Physics, University of Mainz. This work is supported by DFG (SFB443), 
GSI, and the Research Center EMG funded by Forschungsinitiative Rheinland-Pfalz. 

References 

[1] M. Luscher, Comput. Phys. Commun. 165 (2005) 199, hep-lat/0409106; JHEP 12 (2007) Oil, 
arXiv:0710.5417. 

[2] L. Del Debbio, H. Panagopoulos and E. Vicari, JHEP 08 (2002) 044, hep-th/0204125; D.J. Antonio 
et al., Phys. Rev. D77 (2008) 014509, arXiv:0705.2340; S. Schaefer, R. Sommer and F. Virotta, PoS 
LAT2009 032, arXiv:091 0.1465. A. Bazavov et al., Phys. Rev. D81 (2010) 1 14501, arXiv: 1003.5695. 

[3] M. Luscher, Commun. Math. Phys. 293 (2010) 899, arXiv:0907.5491; JHEP 08 (2010) 071, 
arXiv: 1006.4518. 

[4] S. Schaefer, R. Sommer and F. Virotta, Critical slowing down and error analysis in lattice QCD 
simulations, arXiv: 1009.5228. 

[5] M. Liischer, Topology, the Wilson flow and the HMC algorithm, arXiv: 1009.5877. 

[6] S. Capitani et al, PoS LAT2009 095, arXiv:0910.5578. 

[7] B. Knippschild, these proceedings. 

[8] B.C. Tiburzi and A. Walker-Loud, Nucl. Phys. A764 (2006) 274, hep-lat/0501018. 

[9] UKQCD Collaboration (CR. Allton et al.) Phys. Rev. D47 (1993) 5128, hep-lat/9303009. 
[10] A. Hasenfratz and F. Knechtli, Phys. Rev. D64 (2001) 034504, hep-lat/0 103029. 
[11] L. Del Debbio et al, JHEP 02 (2007) 082, hep-lat/0701009. 
[12] L. Del Debbio et al, JHEP 02 (2007) 056, hep-lat/0610059. 
[13] R. Sommer, Nucl. Phys. B411 (1994) 839, hep-lat/93 10022. 

[14] PF Bedaque, Phys. Lett. B593 (2004) 82, nucl-th/040205 1 ; G.M. de Divitiis, R. Petronzio and 
N. Tantalo, Phys. Lett. B595 (2004) 408, hep-lat/0405002; CT. Sachrajda and G. Villadoro, Phys. 
Lett. B609 (2005) 73, hep-lat/041 1033; J.M. Flynn, A. Juttner and CT. Sachrajda, Phys. Lett. B632 
(2006) 313, hepdat/0506016. 

[15] PA. Boyle et al, JHEP 05 (2007) 016, hep-lat/0703005. 

[16] PA. Boyle et al, JHEP 07 (2008) 1 12, arXiv:0804.3971. 

[17] R. Frezzotti, V. Lubicz and S. Simula, Phys. Rev. D79 (2009) 074506, arXiv:08 12.4042. 

[18] F. Jegerlehner and A. Nyffeler, Phys. Rept. 477 (2009) 1, arXiv:0902.3360. 

[19] E. de Rafael, Phys. Lett. B322 (1994) 239, hep-ph/9311316. 

[20] T. Blum, Phys. Rev. Lett. 91 (2003) 052001, hepdat/0212018. 

[21] A. Juttner and M. Delia Morte, PoS LAT2009 143, arXiv:0910.3755; Quark disconnected diagrams 
in chiral perturbation theory, arXiv: 1009.3783. 

[22] B. Jager, A. Juttner, M. Delia Morte and H. Wittig, in preparation. 



7 



